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ABSTRACT 

We use SDO/HMI and SOLIS/VSM photospheric magnetic field measurements to model the force-free coro- 
nal field above a solar active region, assuming magnetic forces to dominate. We take measurement uncertain- 
ties caused by, e. g., noise and the particular inversion technique into account. After searching for the optimum 
modeling parameters for the particular data sets, we compare the resulting nonlinear force-free model fields. 
We show the degree of agreement of the coronal field reconstructions from the different data sources by com- 
paring the relative free energy content, the vertical distribution of the magnetic pressure and the vertically 
integrated current density. Though the longitudinal and transverse magnetic flux measured by the VSM and 
HMI is clearly different, we find considerable similarities in the modeled fields. This indicates the robustness 
of the algorithm we use to calculate the nonlinear force-free fields against differences and deficiencies of the 
photospheric vector maps used as an input. We also depict how much the absolute values of the total force-free, 
virial and the free magnetic energy differ and how the orientation of the longitudinal and transverse components 
of the HMI- and VSM-based model volumes compare to each other. 

Subject headings: Sun: magnetic fields — Sun: photosphere — Sun: photosphere magnetism — Sun: corona 



1. INTRODUCTION 

Presently, no measurements of the magnetic field vector 
in the outer solar atmosphere, the corona, are routinely per- 
formed. Nevertheless, its magnetic structuring is of great in- 
terest and an active field of research. To overcome the unavail- 
ability of coronal magnetic field measurements, photospheric 
magnetic field vector data can be used to reconstruct the coro- 
nal field. Photospheric vector data are obtained from mea- 
surements of Stokes profiles on a regular basis. Based on the 
photospheric field as an input, a force-free field approach, as- 
suming magnetic forces to dominate, can be used to model the 
magnetic field configuration in the atmospheric layers above. 
Those models are characterized by the equations 



V xB = aB 
V -B = 



(1) 
(2) 



and B = B \, s on the photosphere. Here, a is a scalar torsion 
parameter, being constant along field lines because of Equa- 
tion (12). Since, by assumption, the magnetic energy density 
is dominating in the corona, force-free approximations are a 
valid representation of the coronal magnetic field (except for 
the very moment of explosive events). 

Simpler force-free approaches (potential and linear force - 
free, assuming a — and a = constant, respectively) can de- 
scribe the coronal magnetic field topology only to some extent 
since they adhere to restrictive assumptions about the electric 
current. Such models require only the line-of-sight compo- 
nent of the photospheric magnetic field as an input. Large- 
scale electric currents, however, are induced by the reconfig- 
urations of the transver se magnetic field comp onent in which 
excess energy is stored dPriest & Forbesl2002l) . Consequently, 
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potential and linear force-free models can not approximate the 
magnetic energy content of the corona adequately and non- 
linear force-free (NLFF) models are required. NLFF models 
assume a = a(r), i. e., the torsion parameter to be a function 
of space and varying from field line to field line. These are 
currently the most realistic class of force-free fields used to 
model the coronal magnetic field. In particular, it is the ex- 
cess of magnetic energy of a NLFF field over the correspond- 
ing linear force-free field (or a fraction of it), representing an 
upper limit of the "free" magnetic energy that can be either 
transformed or released during eruptions. 

Several studies have focused on the NLFF field reconstruc- 
tion of the coronal magnetic field and its associated magnetic 
energy above solar active regions prior to and after flares. 
Smalle r flares were found to require free energies cc 10 23 J - 
10 24 J (iGilchrist et al.ll20H Uing et al.ll20ll iThalmann et al.1 



120081: IRegnier & Priest! 120071: iRegnier et al.l 120021) while 
large flares seem to occur only when energies ocl0 25 J - 
10 26 J are available (ISun et al.1 120 1 2t Ijmg et all 120 Id 120091 
2008t IThalmann & Wiegelmannll2008t ISchriiver etafl boOoj 
MetcalfetalJl2005l) . 



Force-free field models are applicable for magnetically 
dominated atmospheric layers which, following lGaryl (120011) . 
is true between about 400 km and 200 Mm above photo- 
spheric levels. Thus, vector maps derived from routine mea- 
surements of the Stokes profiles at photospheric levels, repre- 
sent an inconsistent lower boundary condition for the NLFF 
extrapolation. "Preprocessing" of photospheric magnetic field 
vector data is a numerical procedure that is advantageous 
in order to obtain suitable lower bo undary conditions for 
successful force-fre e reconstructions dFuhrmannet al.l 120 lit 
Metcalf et al]|2008l: IWiegelmann et alj|2008t iFuhrmann et al.l 



2007ir Moreover, uncertainties in the measurements, caused, 
e. g., by noise or the particular inversion and azimuth am- 
biguity removal tec hniques, should be taken into account 
(iDeRosa et al.ll200"9l) . The effect of measurement and anal- 
ysis errors on the outcome of NLFF field extrapolations was 
studied by Wie gelmann et alj ( 120101) . They found that mainly 
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the transverse photo spheric field is influenced by instrumen- 
tal effects and that these errors lead to strongest deviations in 
resulting NLFF fields at low heights in the model volume. 

In this work, we use SDO/HMI and SOLIS/VSM pho- 
tospheric magnetic field measurements to model the NLFF 
coronal field above a solar active region. We compare quanti- 
ties like the total magnetic energy content and free magnetic 
energy, the longitudinal distribution of the magnetic pressure, 
and orientation of the magnetic field components in the HMI- 
and VSM-based model volumes. We relate the appearing dif- 
ferences to the photospheric quantities such as the magnetic 
fluxes and electric currents but also show the extent of agree- 
ment of NLFF extrapolations from different data sources. 

2. DATA SOURCES AND EVENT SELECTION 

2.1. Solar Dynamics Observatory - Helioseismic and 
Magnetic Imager 

The Helioseismic and Magnetic Ima ger (HMI) is part o f 
the Solar Dynamics Observatory (SDO; iSchou et all (12012b ) 
and observes the full Sun at six wavelengths and full Stokes 
profile in the Fe i 617.3 nm spectral line. Photospheric line- 
of-sight LOS and vector magnetograms are retrieved from 
filtergrams with a plate scale of 0.5 arc-second (Hoeksema 
T. et al., 2012, in preparation). From filtergrams averaged 
over about ten minutes, Stokes parameters are derived and in- 
ve rted using the Milne-Eddington (ME) inversion algorithm 
of iBorrero et all (120101 ) (the filling factor is held at unity). 
Within automatically iden tified regions of strong magnetic 
fluxes dTurmon et al.l2010l) . the 180° azimuth ambigu ity is re- 
solved using the Minimum Energy Algorith m (MEA; lMetcalfl 
[1991 iMetcalf et al.ll2006t iLeka et alil2009l) . The noise level 
is w 10 G and » 100 G for the longitudinal and transverse mag- 
netic field, respectively. 

2.2. Synoptic Optical Long-term Investigations of the Sun - 
Vector-SpectroMagnetograph 

The Vector-SpectroMagnetograph (VSM) is a part of the 
Synoptic Optical Lon g -term Investigations of the Sun project 
(SOLIS; iKeller et all J2003I) ). The VSM provides, among 
other measurements, full disk vector observations and area 
scans of selected active regions on a daily basis, observing 
conditions permitting. For a typical observing day, full-disk 
photospheric vector magnetograms are available as a quick- 
look (QL) and a ME inversion product dHennev et al.ll2006h 
based on observations in the Fe i 603.15 nm and 603.25 nm 
spectral lines. 

The QL parameters include estimates of the magnetic flux 
density, inclination and azimuth. Its main purpose is to serve 
as an initial guess for the ME inversions. The QL analysis 
estimat es the azimuth and inclination following lAuer et"ai1 
dl977al) . The magnetic flux density is approximated us- 
ing the estimated in clination with the estimated LOS flux 
dRees & Seme il [19791) and assuming that the filling factor is 
one. 

The vector field parameters, for pixels with polariza- 
tion signals clearly above the noise levels, are deter- 
mined using the ME inver sion techniqu e as implemented by 
iSkumanich & Litesl ( 119871) and based on lAuer etaTI dl977bl) . 
The 180° azimuth ambiguity in the transverse field is resolved 
on the entire solar disk using the Non-Potential ma gnetic 
Field Calculation method of (NPFC; iGeoreoulislfcOOSl) . The 
vector data are provided with a plate scale of one arc-second. 
The lower limits for the noise levels are a few Gauss in the 
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Fig. 1.— HMI full-disk LOS magnetic field on 2011 April 12 at 16:46:19 
UT. White/black represents positive/negative magnetic polarities. Yellow 
dashed outlined is the area for which the 180° azimuth ambiguity is resolved 
for the HMI vector map. The blue dotted outline marks the VSM-QL mea- 
surement and the red solid rectangle shows the selected FOV for analysis. 
The color bar indicates the longitudinal magnetic flux density. 

longitudinal and 70 G in the transverse field measurements. 

2.3. Event Selection 

To compare the characteristics of a solar active region, full- 
disk LOS magnetic field observations of the SDO/HMI and 
the SOLIS/VSM were scanned. On 2011 April 12, four ac- 
tive regions, AR 11185, AR 11186, AR 11187 and AR1 1190 
were observed on the solar disk. The latter two, AR 11187 
and AR 11190 appear to be well isolated. For AR 11190 
both HMI and VSM (QL and ME) data are available (see Fig- 
ure [TJ. Vector magnetic field measurements were available 
from SDO/HMI at 16:48 UT and from SOLIS/VSM at 16:43 
UT. We define our field-of-view (FOV) as -404" < x < - 125" 
and 214"<y<381" relative to Sun center, i.e., covering 
about 200 x 120 Mm 2 (red outline around AR 11190 in Fig- 
ure [TJ. The selected FOV is centered around N15E15, with 
the easternmost regions being located at distances not farther 
than about 25° away from disk center. We assume, for the 
purpose of this work, that the vertical and longitudinal field 
are almost identical and that the same is true for the horizon- 
tal and transverse field. Furthermore, we neglect the curvature 
of the Sun and assume a planar geometry of the photospheric 
boundary. 

The purpose of this work is to compare the analysis of the 
photospheric magnetic field and subsequent force-free mod- 
eling based on HMI and VSM (ME and QL) vector maps. We 
investigate if using the not fully ME inverted VSM-QL data 
product significantly changes the quality of the extrapolation. 

3. RESULTS 

3.1. Magnetic structure of the Photospheric Field 

Given the twice as large plate scale of VSM, we bin the 
HMI vector maps to the resolution of VSM in order to com- 
pare the photospheric magnetic field and subsequent force- 
free modeling. We note here that this only slightly affects 
the findings described in the following. The magnetic field in 
AR 1 1 190 has a dipole structure (Figure [2^-c). Similar struc- 
tures are seen in the HMI and VSM data, e. g., the relatively 
compact appearance of the positive and the more fragmented 
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Fig. 2. — Magnetic vector maps from (a) HMI and (b) VSM-ME and (c) VSM-QL. The longitudinal magnetic flux is shown as color-coded background. Blue 
(red) arrows indicate the transverse magnetic field, originating from (positive) negative polarity for regions where B;> 100 G. Contours are drawn at ±200G. 
Histograms of (d) the absolute longitudinal flux fi/, (e) the absolute transverse field B, and (f) the longitudinal current density Ji (calculated at only where 
B, > 100 G and Bi £Q). HMI, VSM-ME and VSM-QL are represented by blue, green and red lines, respectively. Indicated by the dark and light gray shaded 
areas in panel (d) are the "weak" and "strong" field regions for which the longitudinal fluxes are listed in Table[T] 



TABLE 1 

Unsigned longitudinal magnetic flux |0 m | calculated from HMI and VSM 

MAGNETIC VECTOR MAPS. INDICATED ARE THE VALUES FOR TAKING ALL PIXELS INTO 
ACCOUNT ("TOTAL"), CONSIDERING ONLY STRONG FIELD REGIONS (WHERE 

B) > 500 G; "strong") and only weak field regions (where 
100 G < Bi < 500 G; "weak"). Also given are the values of the net Lorentz 

FORCE (fforce) AND NET TORQUE (^torque)- 



Data 


10m 


|[X10 22 


Mx] 


^force 


^torque 




total 


strong 


weak 




HMI 


1.46 


0.57 


0.62 


0.49 


0.44 


VSM-ME 


2.03 


0.65 


1.09 


0.82 


0.87 


VSM-QL 


1.35 


0.17 


0.81 


0.78 


0.85 



pattern of the negative polarity region. The transverse field, 
B,, has a clearly diverging (converging) structure above the 
positive (negative) polarity. 

The unsigned magnetic flux, when calculated from 
the HMI (|< MI |), VSM-ME (|^ SM " ME |) and VSM-QL 
(l0m SM ^ L |) data is balanced to within ss 10%. All are 
on the order of 10 22 Mx, where |^ MI | ~ 0.7 |<7^ SM ~ ME | and 
|< MI | -1.1 |<7)m SM ~ QL | (see Table [B- The net flux of the HMI 
and VSM vector maps is positive, (p m cc 10 21 Mx and about 
10% of the unsigned flux. In the VSM-ME vector maps, zeros 
are assigned to pixels where the measured polarization signal 
is too weak to perform a reliable inversion. If we disregard 
these "zero-pixels", about 60% of the total number of pix- 
els in the HMI and VSM vector maps remain for compari- 
son. Then the unsigned flux in the HMI (VSM-QL) vector 
map is a; 3% less and the net longitudinal flux is » 1% (« 4%) 
less (more). In other words, the unsigned and net longitudinal 
fluxes of the HMI and VSM-QL vector maps are practically 
unchanged when disregarding the zero pixels in the VSM-ME 
vector map and their contribution is negligible. 

In the following, we study the spatial distribution of the flux 
and consider regions of strong and weak magnetic flux sepa- 
rately, defining "strong" as where B\ greater than 500 G and 
"weak" as where 100 G < 5/ < 500 G (where B, = |Blos|). The 



zero-pixels do not cover parts of the weak and strong field re- 
gions, i. e., they don't influence the calculated quantities. For 
all vector maps the weak flux regions have more flux when 
compared to the strong flux regions: « 10% more in HMI, 
* 40% more in VSM-ME and * 80% more in VSM-QL (the 
latter being partly a consequence of the cutoff of ^vsm-ql ^ 
ss 1 kG). Weak field regions cover 10% and 20% of the FOV 
of HMI and VSM, respectively. In these regions the VSM 
data have more longitudinal flux (|^ SM " ME | ~ 1.7 |< MI | and 

\<Pm M ^ L | ~ 1 -3 |0™'D- The strong field regions cover around 
3% of the considered area of the HMI and VSM-ME vec- 
tor maps, and 1% of the VSM-QL data. In those regions, 
the HMI (VSM-QL) data reveal only about 90% (25%) of 
|^vsm-mE| Com p ar ed to HMI, VSM-QL has less longitudi- 
nal flux in strong field regions but more in weak field regions. 
This might partly be attributable to the influence of scattered 
light. The relative portion of scattered light is higher in strong 
field regions so that the effect on the measured flux is larger, 
i. e., measured fluxes are more strongly reduced. This level 
of scattered light in the VSM vector maps is less than 10%, 
while its importance for HMI vector maps is not quantified 
yet. 

The distribution of Bi decays rapidly with increasing flux 
density (see Figure |2}l). On overall, the longitudinal field of 
the HMI and VSM-ME vector maps are comparable, espe- 
cially at higher field strengths. In general, VSM-ME detects 
more longitudinal field than HMI. While b^sm-me does nQt 

exceed ~ 1.8 kG, Bf m can be as high as * 2.1 kG. A differ- 
ent behavior is found for the VSM-QL data. The HMI and 
VSM-QL estimates match at field strengths below as 0.7 kG, 
the relative occurrence of ^vsm-ql clearly differs at higher 
values. A similar trend is found when considering the abso- 
lute value of the transverse field, B,-\B,\, (Figure^). Maxi- 
mal values of B^sm-me are below « 1.3 kG while B™ 1 can be 
up to « 1.8 kG. The shapes of the B, distributions differ with 
larger discrepancies at higher values of B,. Only at values 
0.fkG<B,<0.8kG the distributions of £ v sm-mi? and b hmi 
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Fig. 3. — Preprocessed magnetic vector maps from (a) HMI and (b) VSM-ME and (c) VSM-QL. The longitudinal magnetic flux is shown as color-coded 
background. Blue (red) arrows indicate the transverse magnetic field, originating from (positive) negative polarity for regions where B/> 100 G. Contours are 
drawn at ± 200 G. Histograms of (d) the absolute longitudinal flux B/, (e) the absolute transverse field B, and (f) the longitudinal current density J; (calculated at 
only where B, > 100 G and B/ # 0). HMI, VSM-ME and VSM-QL are represented by blue, green and red lines, respectively. Indicated by the dark and light gray 
shaded areas in panel (d) are the "weak" and "strong" field regions for which the longitudinal fluxes are listed in Table[2] 



TABLE 2 

Unsigned longitudinal magnetic flux \ip m \ calculated from the 

PREPROCESSED HMI AND VSM MAGNETIC VECTOR MAPS. INDICATED ARE THE 
VALUES FOR TAKTNG ALL PIXELS INTO ACCOUNT ("TOTAL"), CONSIDERTNG ONLY 
STRONG FIELD REGIONS (WHERE B) > 500 G; "STRONG") AND ONLY WEAK FIELD 

regions (where 100 G < £/ < 500 G; "weak"). Also given are the values of 

THE NET LORENTZ FORCE (£f olcc ) AND NET TORQUE (Etorque)- 



Data 


10m 


|[xl0 22 


Mx] 


^force 


^torque 




total 


strong 


weak 




HMI 


1.42 


0.41 


0.71 


0.01 


0.01 


VSM-ME 


2.03 


0.62 


1.11 


0.01 


0.02 


VSM-QL 


1.35 


0.16 


0.81 


0.01 


0.02 



are relatively close to each other. Similar to the trend of the 
longitudinal field, the transverse field estimates of VSM-QL 
increasingly differ with increasing B, and, on overall, VSM- 
QL detects less transverse field. 

We also show the absolute value of the longitudinal current 
density, 7/ = |J/|, distribution in Figure |2f. We calculate the 
longitudinal current density J; everywhere where B, > 100G 
and B) + 0. The // distribution clearly differ over the entire 
range, most pronounced at lOmAm <//< 30mAirT 2 . On 
overall, J\ is strongest in the HMI vector maps and weakest 
in the VSM-QL data. The mean values of J\ are 5 mAirT 2 , 
0.6 mArrT 2 and 0.7 mArrT 2 , respectively. The spectral power 
of // is more or less constant on all spatial scales in the HMI 
vector maps. For the VSM vector maps more spectral power 
of /; is found on larger scales than on small scales. The total 
longitudinal current is approximately -4 x 10 7 A, -2 x 10 7 A 
and -1 x 10 7 A for the HMI, VSM-ME and VSM-QL vector 
maps, respectively. 

We also probe the influence of the different methods used 
to resolve the 180° azimuth ambiguity in the transverse field 
of the VSM and HMI vector maps. We compare B/, B, and J\ 
of the MEA ambiguity resolved HMI vector maps to that of 
NPFC ambiguity resolved HMI maps (courtesy of M. Geor- 
goulis). The power spectra of B/, B, or J\ as a function of 



wave number (i. e., length scale) are very similar and show no 
significantly different distributions. The power of Bi and B, 
decreases with increasing wave number (i. e., with decreas- 
ing length scale, prominently on scales < 20 Mm) while the 
power of J i if well distributed over all length scales. We there- 
fore assume that the method used to resolve the 180° azimuth 
ambiguity has a minor impact on B[, B, and 

All of the above indicates that most longitudinal magnetic 
flux is detected by VSM-ME, while the transverse field, and 
hence electric currents, are partly of spatial scales which 
VSM-ME cannot fully resolve. Weakest signals for all, B/, 
B, and //, is found in the VSM-QL data. However, as will 
be shown below, the relative occurrences of all quantitites are 
clearly more similar after applying a "preprocessing" routine. 

3.1.1. The Effect of Preprocessing 

Vector maps, suitable as input for our force-free model 
technique are expected to be flux-balanced with a van- 
ishing net force and net torque. We use dimensionless 
parameters, ef orce and etomue, to quantify such properties 
dWiegelmann & Inhesteril2006l:lATvlll989tlMolodenskvll 19741: 
iMolodenskiil 1 1 9691) . For AR 11190, ef orC e and etorque are on 
the order of 1.0 (see Table[TJ. These are similar to the values 
calculated from vector data used in previous studies (Wiegel- 
mann T, et al., 2012, Sol. Phys., submitted; Thalma nn et alj 
120081) . showing that the measured data are, in principle, 
inconsistent with the force- free field approach. We apply 
the preprocessing routine of Wie gelmann & Inhesterl (120061) 
which reduces the net force and net torque. Doing this, we 
do not intend to suppress the existing forces in the photo- 
sphere (where we know that the field is not force-free). In- 
stead, we try to approximate the magnetic field at a chromo- 
spheric level where magnetic forces are expected to be much 
smaller than in the layers below. From the preprocessed vec- 
tor maps, which we assume to be situated somew here betwee n 
400 km and 800 km a bove the photosphere (e. g jGarvl (120011) : 
iMetcalfetal] ( fl995l) ). we find 6f orce oc 10 2 and e torque oc 10 2 
(see Table©. 
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The preprocessing influences the structure of the magnetic 
vector data. As can be seen in Figure |3^-c, the prepro- 
cessed B, of HMI, VSM-ME and VSM-QL are more simi- 
lar than before applying the preprocessing (as shown in Fig- 
ure|2^-c). This is because the preprocessing scheme, not only 
smooths B, but also alters its values in order to reduce the 
net force and torque. This is physically motivated since B, is 
measured with lower accuracy than the longitudinal magnetic 
field. The effect of smoothing can be recognized by compar- 
ing the ± 0.2 kG contour of the original and preprocessed HMI 
vector map (Figure |2K and Figure (3Jl, respectively). It is less 
fragmented in the preprocessed data, compared to that in the 
original data. 

As in the original data, the preprocessed HMI, VSM-ME 
and VSM-QL data have more longitudinal flux in the weak 
field regions than in the strong fields (»30%, «40% and 
« 80% more, respectively; see Table As for the raw 
data, the strong field regions cover about 2%, 3% and 1% 
of the considered area of HMI, VSM-ME and VSM-QL. 
In those regions, the HMI data reveal ^80% and VSM- 
QL «30% of \4>l SM - ME \. Weak field regions cover about 
15% and 20% of the FOV of HMI and VSM, respectively. 
In these regions, the VSM-ME data have most longitudinal 
flux (|< SM ~ ME I ~ 1.6 |< MI |), followed by the VSM-QL data 

(|^ n SM - QL |-1.2|< MI |). 

The relative occurrences of Bi in the preprocessed HMI, 
VSM-ME and VSM-QL vector maps are clearly more consis- 
tent with each other than for the raw data (Figure |3jl)- This is 
to be expected since the preprocessing scheme only smooths 
the longitudinal field (i. e., does not account for any flux im- 
balance) while it smooths and alters the transverse field. The 
latter obviously results in the histograms of B, agreeing well 
over the entire range of strengths after the preprocessing ( Fig- 
ure^)- We suspect that this is due to smoothing of the in- 
strumental noise. The mean value of the changes due to pre- 
processing in the longitudinal field is oclO~ 7 G and for the 
transverse field on the order of 10 G (about 1% of the peak 
transverse flux). Therefore, the mean changes due to the pre- 
processing are well within the measurement uncertainty of the 
HMI and VSM (though the largest local modifications can be, 
for both B\ and B,, several 100 G). The modification due to 
preprocessing also includes the rotation of B, (i. e., both B x 
and B y change sign). This affects 2.0%, 3.2% and 0.8% of the 
locations where B, > 250 G in the HMI, VSM-ME and VSM- 
QL vector maps, respectively. For pixels with transverse fields 
on the order of the measurement error or less, the field is ro- 
tated in 6.9% of the cases in the HMI and VSM-ME vector 
maps and in 3.3% of the cases in the VSM-QL data. By coin- 
cidence, this rotation could cause a flipping of the transverse 
field (i. e., a rotation of 180°) but this was not observed when 
comparing the original with the preprocessed data. 

Altering the transverse magnetic field components by pre- 
processing necessarily leads to a change in the electric cur- 
rents. The relative occurrence of Ji (Figure |3f) shows that the 
currents of the preprocessed HMI and VSM-ME vector maps 
are very similar. Comparing, Figures 01 and [3}l, it can be 
seen that the highest absolute values // decrease, about as 1 1% 
for HMI and about * 18% for the VSM-ME. The mean val- 
ues of J/ are «0.9mArrT 2 , *0.6mAm -2 and =s0.4mAirT 2 , 
respectively. Though reducing the total longitudinal current, 
the preprocessing does not redistribute the current to larger 
scales, it causes a decrease of the spectral power of 7/ at small 
scales (especially on scales below » 6 Mm). The total longi- 



tudinal current decreases to about -1 x 10 7 A, -3 x 10 6 A and 
-1 x 10 6 A in the preprocessed HMI, VSM-ME and VSM-QL 
vector maps, respectively. 

3.2. Force-free Modelling 

As pointed out by De Rosa et al.l ( 120091) . even after prepro- 
cessing the photospheric vector maps, force-free model tech- 
niques might still not find a strictly force- and divergence-free 
solution in the coronal volume above. Therefore, we calcu- 
late the NLFF fields using an e xtended optimization scheme 
dWiegelmann & Inhesterfl20lol) . By construction, the NLFF 
field drops towards a potential field, prescribed on the lat- 
eral and top boundaries of the volume-bounding surface, &V, 
within a finite-size layer during the iterative solution pro- 
cess. For further computations and subsequent analysis, we 
extract the inner physical domain (that is the sub-volume, ex- 
cluding the boundary layer towards the side and top bound- 
aries). The resulting domain of analysis is of the dimension 
'Vsub ~ 179 x 98 x 75 Mm 3 . We use the normal component of 
the magnetic field, B„, on the volume-bounding surface of this 
NLFF sub-field, d r V su h, to calculate the corresponding poten- 
tial field Bp 0t = Vi^, where iff is the solution to the Laplace 
equation, subject to d„if/ = £>" lff on d'V^. To quantify the 
quality of the reconstructed NLFF fields, we list the normal- 
ized mean Lorentz force (-CforceX the mean field divergence 
(-Cdiv) an d the curren t- wei ghted mean angle, 9j, between B 
and J (see Appendix lA.il ) in Table [3] 

The improved optimization scheme allows us to relax the 
magnetic field also on the lower boundary (besides minimiz- 
ing the volume-integrated force and divergence; see Appendix 
|AV The relaxation of the lower boundary introduces a fur- 
ther modification of the vector data, in addition to that by the 
preprocessing applied before. The mean value of the change 
of the longitudinal fields due to the relaxation of the lower 
boundary is oc 10~ 4 G and largest local modifications are of 
the order 1 G. The mean changes of the transverse field are on 
the order of 10 G and largest local modifications can be sev- 
eral 100 G. Given the noise levels of HMI and VSM measure- 
ments of the longitudinal (w 10 G and a fewG, respectively) 
and transverse field (w 100G and >70G, respectively), the 
mean (largest local) modifications are well within (on the or- 
der of) the measurement error. In addition to the modification 
of the absolute fluxes, the transverse field vector is rotated in 
some places (i. e., the sign of both B x and B y changes). For 
transverse fields stronger than 250 G, such a rotation occurs 
in only about 0.1% of the corresponding pixels in the pre- 
processed VSM data but nowhere in the preprocessed HMI 
vector map. The rotation predominately occurs in regions of 
weak transverse fields close to the measurement uncertain- 
ties (at 17.7%, 14.7% and 1.7% of the pixels with B,< 100 G 
in the preprocessed HMI, VSM-ME and VSM-QL data, re- 
spectively). Comparing the preprocessed vector maps and 
the lower boundary after optimization, no places were found 
where a flipping (i. e., a rotation by 180°) occurred. 

As can be seen in Table[3](see also Appendix lA.2l ), the en- 
ergy of the potential and the NLFF field, extrapolated from 
the VSM-ME data is a factor of a 2 higher than that of the 
model based on HMI or VSM-QL data. This is also true for 
the upper limit for the free magnetic energy, AE m . The un- 
signed longitudinal flux on the VSM-ME bottom boundary 
has * 30% more unsigned longitudinal flux than the HMI and 
VSM-QL lower boundary. Since the longitudinal field com- 
ponent enters in form of boundary conditions in the poten- 
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TABLE 3 

VlRIAL (£ v [ r ) AND TOTAL MAGNETIC ENERGY (E n [ff ) OF THE NLFF FIELD, POTENTIAL 
FIELD ENERGY (£pot) AND THE UPPER LIMIT FOR THE FREE MAGNETIC ENERGY 

(A£" m = £„iff - Spot). Normalized mean Lorentz force (Xforce), mean field 

DIVERGENCE (Xdiv) AND CURRENT-WEIGHTED MEAN ANGLE, 6j. 



Data 


.Evil - 


E n \s £pot 
[x 10 25 J] 


A£ ra 


-£force 


-£div 


0j 

n 


HMI 


2.97 


2.68 2.49 


0.19 


0.43 


0.24 


3.8 


VSM-ME 


4.12 


5.14 4.55 


0.59 


0.54 


0.33 


5.2 


VSM-QL 


1.89 


2.14 1.85 


0.29 


0.50 


0.33 


4.4 



tial field calculation, it is not surprising that the VSM-ME 
field model has more potential energy. The energy of the 
NLFF field is expected to be higher than the potential en- 
ergy, since it also incorporates the transverse field compo- 
nents. This is true for all models. The relative amount of 
energy, available to be set free during reconfigurations of the 
magnetic field, AE m , is «7%, « 11% and * 14% of the re- 
spective NLFF field energy, E^s, for the HMI, VSM-ME and 
VSM-QL models. Based on the preprocessed lower bound- 
aries, we also calculate the to tal force-free energy g iven by the 
magnetic virial theorem (see iMetcalf et al.l (12005)) and Equa- 
tion |A6| in Appendix |A.2| i. We find a virial energy oc 10 25 J 

(£ VSM-ME >jB HMI >jB VSM-QL ) ^ comparable tQ ^ In 

particular, we find that for the HMI (VSM) models, the virial 
energy is higher (lower) than the NLFF field energy. 

The absolute value of the longitudinal electric current den- 
sities, Jj, integrated over the lowest 36 Mm in the NLFF field 
models based on HMI, VSM-ME and VSM-QL data is shown 
in Figure 2^-c. (We calculate the current density only at lo- 
cations in the volume where B, > 100 G and Bi ^0.) Their 
origin in the underlying positive (20 Mm < x < 70 Mm) and 
negative (90Mm<x< 130 Mm) magnetic polarity regions is 
clearly visible. The HMI-based integrated J\ shows more fine 
structures than that based on VSM data. Above the positive 
polarity area we find a, more or less, radially outward directed 
current structure. Above the negative polarity region, the 
current-related structures are more diffusely distributed which 
might be due to the underlying photospheric field being more 
fragmented. There are also strong, isolated current patches 
in the VSM-ME model (e. g., in the north-eastern part of the 
AR, above the positive polarity region) which are weaker in 
the HMI and VSM-QL models. On overall, the pattern of the 
longitudinally integrated current density is the same for the 
HMI and VSM model fields. 

We also investigate the magnetic pressure, p m , in a longitu- 
dinal cross-section at y = 53 Mm (as indicated by a horizontal 
line in Figure |4^-c) of the NLFF fields (Figure EJl-f). The 
overall pattern of p m appears to be the same when calculated 
from the HMI (p™ 1 ), VSM-ME (/^ SM " ME ) and VSM-QL 

(Pm SM ^ L ) NLFF model volume. Again, correspondence to 
the underlying positive and negative polarity regions is vis- 
ible. On overall, /^ SM - ME > p™ 1 > pl SM ~ QL for same loca- 
tions in the cross section. This is expected, since p m is pro- 
portional to E m , for which a similar trend was found. Iso- 
contours, indicating the magnitude of /?^ SM ~ ME , are located 

higher up in the atmosphere compared to p™ 1 and p^ M ~^ L , 
indicating that the VSM-ME NLFF field is sustained by a 
higher magnetic pressure from below. Three sub-systems are 
clearly identifiable above the positive polarity region, around 
x * 50 Mm, most pronounced in p™ 1 and subsequently more 

diffuse in /v SM ^ and p^ SM ~ ME . A longitudinal dip-like 



structure is clearly visible in p^ SM ME and p^ M ^ L at al- 
titudes z> 10 Mm but very weak in p™ 1 . Above the neg- 
ative polarity region, another two-component sub-system is 
revealed, above 105 Mm < x < 140 Mm. 

To visualize the alignment of the magnetic field, B, and the 
electric current, J, we show cos(^) = (J • B)/(B J) in the lon- 
gitudinal cross-section at y — 53 Mm (Figure |4g-i). A value 
of zero means that the current is perpendicular to the mag- 
netic field. A value of (-)+l indicates that the current is 
(anti)parallel to the magnetic field, as it should be for a per- 
fectly force-free field. On overall the same basic alignment 
pattern of cos(^) is found in all model fields though they are 
shifted and deformed when compared to each other. Some 
distinct features are clearly recovered, e.g., a low-lying arc- 
like structure (region "a" in Figure |4g-i). Similarly, a dip- 
like structure as well as a vertically oriented pattern are dis- 
cernible in all models (regions "b" and "c" in Figure |4g-i, 
respectively). However, there are also places where the elec- 
tric current is parallel to the HMI field while it is anti-parallel 
to the VSM-ME and/or VSM-QL field. Those places can- 
not be attributed to locations of particularly weak currents 
and we suppose that these are consequences of the differing 
non-potential (current-carrying) contributions to the NLFF 
field solutions (see below). The appearing discrepancies may 
partly be caused by not looking at exactly the same topologi- 
cal structure in the different models when displaying at a cer- 
tain longitudinal cross-section of the computational volume. 

We investigate the similarity of the VSM-ME, VSM-QL 
and HMI mod els by comparing the vector correlation (see 
Appendix |A.2t of the potential fields (C pot ), the NLFF fields 
(Cnifr) as well as of the non-potential part (current-carrying; 
B curr = B n is - Bp t), C cun -. In all cases, C pot > C n ifr > C curr 
with Cp t - 1 (i.e., the potential fields are almost identical) 
and Cniff > 0.95 (see Figure [5J. The average vector correla- 
tion between the potential fields based on the HMI and VSM- 
ME data as well as between the HMI and VSM-QL data is 
0.98 (dotted lines in Figure [5^ and |5j), respectively). The 
potential fields based on the VSM-ME and VSM-QL data 
agree best with C pot - 1 (dotted line in Figure [5};). The aver- 
age vector correlation between the NLFF fields is 0.97, 0.94 
and 0.98 for HMI vs. VSM-ME, HMI vs. VSM-QL and VSM- 
ME vs. VSM-QL, respectively. The high values of C pot and 
C n iff are not surprising since the NLFF fields do not devi- 
ate too much from the corresponding potential fields. This 
can be seen from the relative amount of free magnetic energy 
of only as 10% for all models. The vector correlation of the 
current-carrying parts only, however, shows a different behav- 
ior. On average, the vector correlation between the current - 
carrying parts is 0.53, 0.32 and 0.66 for HMI vs. VSM-ME, 
HMI vs. VSM-QL and VSM-ME vs. VSM-QL, respectively. 
A varying pattern with height is found, with a poor correla- 
tion at low altitudes (below « 40 Mm) and steadily increasing 
for higher altitudes. The increasing values with increasing al- 
titude reflect the nearness of the NLFF fields to the potential 
field solutions. With increasing altitude, the correlation of the 
potential field solutions is more and more dominant since, by 
construction, a NLFF field drops to the potential field towards 
the top boundary. To summarize, though the potential and 
NLFF field models based on the HMI, VSM-ME and VSM- 
QL data are similar and show a high vector correlation, this is 
not true for the current-carrying part of the fields where strong 
deviations, especially at low altitudes, are found. 

4. CONCLUSION 
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Fig. 4. — Vertically integrated current density (where B, > 100 G), computed from the NLFF extrapolated fields of (a) HMI, (b) VSM-ME and (c) VSM-QL, 
over the lowest 36 Mm. Black (gray) contours of the longitudinal field at the bottom of the model volume are drawn at +(-)200 G. Magnetic pressure p m below 
36 Mm in the longitudinal cross-section at y = 53 Mm for the (d) HMI, (e) VSM-ME and (f) VSM-QL NLFF models. Panels (g), (h) and (i) show the cosine of 
the angle between the magnetic field and the electric current in the same longitudinal cross-section. The footprint of the longitudinal cross-section on the lower 
boundary is indicated by the black horizontal line in Figuref4t-c. 
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Fig. 5. — Vector correlation of the model fields as a function of height in the volume of interest. The correlation between the potential fields, Cp 0l , is shown 
as dark gray dotted lines. The correlation of the NLFF fields, C n iff, is drawn as black solid lines. That of the current-carrying part of the fields (C cu rr; with 
B cun - = Sniff ~ Sp t) is represented by the light gray dashed-dotted lines. The comparison HMI and VSM-ME models is shown in (a), that of the HMI and 
VSM-QL models in (b), and that of the VSM-ME and VSM-QL models is displayed in (c). 



We use photospheric vector maps, observed on 201 1 April 
12 to model the nonlinear force-free field above AR 11190. 
We use the ME and QL products from SOLIS/VSM as well 
as SDO/HMI data. We compare the measured photospheric 
VSM-ME, VSM-QL and HMI data. We find that VSM-ME 
detects most longitudinal flux, followed by HMI and VSM- 
QL. In strong and weak field regions VSM-ME detects more 
longitudinal flux than HMI. This is not true for VSM-QL: it 
shows less flux than HMI in weak field regions but more flux 
in strong field regions. HMI is found to detect most trans- 
verse field. The transverse fields are not detected as well in 
the VSM data. On overall, the longitudinal electric current is 
strongest in the HMI data, followed by that of VSM-ME and 
VSM-QL. 

Applying a preprocessing routine to construct boundary 



conditions, compatible with the force-free approach, results 
in the overall structure of the transverse field of the HMI and 
VSM vector maps becoming more similar. The longitudinal 
and transverse magnetic flux of the preprocessed VSM-ME 
and HMI vector maps agree very well as does the longitu- 
dinal electric current. Only the VSM-QL data has a clearly 
lower signal in all photospheric measures. Using the prepro- 
cessed HMI and VSM vector maps as an input, we model the 
force-free field above AR 11190. Thus, we are able to quan- 
tify the effect of the choice of the lower boundary condition 
on the model outcome. 

The different longitudinal magnetic flux translates to dif- 
ferences in the force-free field energies of the modeled coro- 
nal volumes. The NLFF, potential and, consequently, the free 
magnetic energies (all oc 10 25 J) of the VSM-ME model are 
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a factor of » 2 higher than those of the HMI and VSM-QL 
models. The models' free energy is between w 7% and * 14% 
of the NLFF field energy, corresponding to «2 - 3 x 10 24 J. 
The virial energy is comparable to the NLFF field energy and 
slightly (higher) lower in the (HMI-) VSM-based models. 

It is believed that flare sizes relate to the available amount 
of free energy in the coronal volume. Only one larger event 
(a Ml. 3 flare on April 15, peaking at 17:12 UT) originated 
from AR 11190. A smaller C 1.1 flare took place on April 
13 and a CI. 6 flare on April 14. From our force-free models, 
we find a free magnetic energy oc 10 24 J on April 12. This 
is cert a inly enough to power C -class flaring (e. g.. Uing et al.l 
( 120101) ; iThalmann etal] (120081) ). The question arises if this 
amount of energy plus the additional built up energy during 
the next two days were sufficient to power the M-class flare 
about three days later. It cannot be ruled out that this amount 
of pre-flare energy is sufficient to power the event but it seems 
rather low. 

Agreements of the HMI- and VSM-based extrapolations are 
also found, e.g., in the vertical distribution of the magnetic 
pressure and the longitudinally integrated current density. The 
longitudinal and transverse vector components of the recon- 
structed NLFF fields match each other to quite high degree 
when viewed as function of height in the atmosphere. The po- 
tential and NLFF field models based on the HMI, VSM-ME 
and VSM-QL data are similar and show a high vector correla- 
tion. This is not true for the current-carrying (non-potential) 
part of the model fields where strong deviations, especially 



at low altitudes, are found. The orientation of the electric 
current density with respect to the magnetic field in a verti- 
cal cross-section of the extrapolation volume is found to be 
similar in portions of the model volumes. However, also dis- 
crepancies are found for some places where, e. g., the current 
is parallel to the field in the model based on the HMI data 
but anti-parallel in the models based on VSM data. On over- 
all, the co-alignment pattern of the electric current and mag- 
netic field in the different models is distorted which may be 
attributed to the differing non-potential portions of the NLFF 
model volumes. 

We conclude that NLFF field extrapolations based on HMI, 
VSM-ME, VSM-QL vector maps result in comparable coro- 
nal field models and related quantities. Most importantly, the 
models agree well regarding the estimated relative free mag- 
netic energy which can be released during explosive events. 
This allows us to assume that our method to reconstruct the 
coronal magnetic field is robust to deficiencies in the lower 
boundary data. 

SOLIS data are produced cooperatively by NSF/NSO and 
NASA/LWS. SDO data are courtesy of the NASA/SDO HMI 
science team. J. K. T. acknowledges supported by DFG grant 
WI 3211/2-1, X.S. is financed by NASA contract NAS5- 
02139 and T. W. is funded by DLR grant 50 OC 0904. We 
thankfully acknowledge M. Georgoulis for ongoing collabo- 
ration and helpful discussions, and the careful consideration 
of the referee to improve the manuscript. 



APPENDIX 

A. NONLINEAR FORCE-FREE MAGNETIC FIELD MODELING 

The improved optimization scheme of IWiegelmann & Inhesterl (120101) . based on a proposed principle by Wheatland et al.l 
( 2000) and originally implemented and extended by [Wiegelma nnl (120041) solves the force-free equations by minimizing 

£o P t = f (w f B- 2 \(V x B) x B\ 2 + w d |V ■ B\ 2 ) dTV + v jf( B - B obs ) ■ W ■ ( B - B ohs ) dS. (Al) 

where B-\B\ and Wf and are position-dependent weighting functions. Wf and wa are set to 1 in the volume of interest and 
drop to zero in form of a cos-profile within a finite size boundary layer towards the lateral and top boundaries. This reduces 
the influence of the boundaries of the computational cube on the solution inside. The boundary layer spans 16 pixels in the 
extrapolation volumes with the HMI vector maps being binned to the resolution of VSM. 

The major improvement of our optimization method is the introducing of the surface term in (1A1K i. e., to relax the field also 
on the bottom boundary. B and B b s denote the iterated and observed magnetic field on the lower boundary of the extrapolation 
volume. Before relaxation, the obs erved photospheric field is pr eprocessed in order to achieve suitable force-free consistent 
boundary conditions (for details see Wiegel mann & Inh ester 2006). The free parameters controlling the net force and net torque 
are set to jUi = 1 and hi = 1, respectively. The deviation from the original data is weighted with ^ = 10~ 3 and the level of 
smoothing is set to jj.4 = 10~ 2 . 

The speed with which the lower boundary is injected during the extrapolation is controlled by v. The elements of the newly 
introduced space-dependent diagonal error matrix W are defined inversely proportional to the squared measurement error of 
the respective field component. Ideally, one would define its elements based on the error introduced by, e.g., the method of 
inverting the measured Stokes profiles or algorithm used to resolve the 180° azimuth ambiguity of the transverse field. Such error 
estimations are, however, not (yet) available so we have to make reasonable assumptions on the nature of W. Therefore, we chose 
its element controlling the LOS fiel d component, Wlqs, to b e unity and the one controlling the transverse field, Wtians> as small 
but positive. As in a previous study (IWiegelmann et al. 2012), we tested different choices of W and v. In particular, we tested all 
combinations of Wtrans K [B t ,B 2 , y[B~,~\ and v- [ 10 _1 , 10~ 2 , 10~ 3 ], where B,-\B,\. We found the NLFF field, calculated based on 
the choice w tra ns K B 2 and v = 10~ 3 , performs best for both HMI and VSM input data (for the quantification of the performance 
see Section lATb . 

A.l. Quality of the Nonlinear Force-free Magnetic Field Modeling 

Three measures to quantify the goodness of the reconstructed fields are used. These are the mean Lorentz force (Xforce) and 
field divergence (Xdiv) and the current weighted mean angle between the magnetic field and the electric current density (6>j), 
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which are of the form 



_ v \JjX Bi\ 2 

■i-foice ~ / , B 2 
i i 



6 j = sin" 1 



I 



Ui x Bj\ 
Bi 



(A2) 
(A3) 

(A4) 



where / = \ J\. -Cf OIce and Xdiv are normalized to a unit volume. For a perfectly force-free case, all the above quantities vanish. 
For the extrapolations presented in this study, these values are -Cforce = 0.43, -Cdiv = 0.24 and 6j = 3.8° for HMI. For the VSM-ME 
(-QL) NLFF models they are £ force = 0.54 (0.50), -C d iv =0.33 (0.33) and 0j=5.2° (4.4°). All metrics are significantly better 
compared to the corresponding values when not relaxing the lower boundary. For comparison, one then finds -Cfoice =23.5, 
£ d i v = 12.4 and 0j = 23.2 for HMI and £ forc e = 26.2 (26.3), £ d iv = 17.8 (18.8) and 8j = 37.3° (36.8°) for VSM-ME (-QL). 

A.2. Calculation of Related Quantities 
The magnetic energy content of the modeled NLFF fields is calculated as 



— f 

2/io Jv 



B d'V. 



(A5) 



An upper limit for the free magnetic energy can be estimated by subtracting the magnetic energy of a potential field (E pot ) from 
the energy of a NLFF field (E a \s), i- e -> A£ m = E n if{ - E pot . 

The magnetic virial theorem can be used to calculate the force-free field energy in a volume solely from its boundary. When 
applied to the lower boundary, it is given by: 



/4> Js 



(xB x + yB y )B z dS. 



(A6) 



Strictly speaking, one has to take all the magnetic flux at the lower boundary into account, i.e., no field lines are supposed 
to connect to the outside of the NLFF field. Furthermore, E vlI depends on the origin of the coordinate system unless the net 
transverse Lorentz forces vanish (IWheatland & Metcalfll2006l) . 
To quantify the degree of agreem ent between two given vector fields, B and £?', we use the "vector correlation" metric, as 

ment of thi 

C(B,B') = = ' ' . (A7) 



lo quantity the degree of agreement b< 
introduced by lSchrijver et al. (2006). The agreement of the two vector fields is defined as 



where £?, and B. are the two arbitrary vectors fields at each point i. For two identical vector fields £?, = B' r C(B, B') = 1 
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